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^ . Abstract 

We study the evolution of the flavour non-singlet deep-inelastic structure functions i^NS 
and F3 at the next-to-next-to-next-to-leading order (N 3 LO) of massless perturbative QCD. 
The present information on the corresponding three-loop coefficient functions is used to 
derive approximate expressions of these quantities which prove completely sufficient for 
values x > 10~ 2 of the Bjorken variable. The inclusion of the N 3 LO corrections reduces 
the theoretical uncertainty of a s determinations from non-singlet scaling violations arising 
from the truncation of the perturbation series to less than 1%. We also study the pre- 
dictions of the soft-gluon resummation, of renormalization-scheme optimizations by the 
principle of minimal sensitivity (PMS) and the effective charge (ECH) method, and of the 
Pade summation for the structure-function evolution kernels. The PMS, ECH and Pade 
approaches are found to facilitate a reliable estimate of the corrections beyond N 3 LO. 



1 Introduction 



Structure functions in inclusive deep-inelastic lepton-nucleon scattering (DIS) are among 
the observables best suited for precise determinations of the strong coupling constant a s . 
At present their experimental uncertainties result in an error A exp o; s (M§) ~ 0.002 at 
the mass of the Z-boson |L]. A further reduction of this error can be expected, especially 
from measurements at the electron-proton collider HERA after the forthcoming luminosity 
upgrade. The standard next-to-leading order (NLO) approximation of perturbative QCD 
summarized in ref. 0, on the other hand, leads to a theoretical error A t ^a s (M^) « 0.005. 
This error is dominated by the uncertainty due to the truncation of the perturbation series 
as estimated from the renormalization-scale dependence. Hence calculations beyond NLO 
are required to make full use of the present and forthcoming data on structure functions. 

The ingredients necessary for next-to-next-to-leading order (NNLO) analyses of the 
structure functions in Bjorken-x spaceQ have not been completed up to now: Unlike the 
two-loop coefficient functions which were calculated some time ago (and completely 
checked recently ||), only partial results 0, §, [| [10|] have been obtained for the three- loop 
splitting functions so far. However, we have recently demonstrated |L1[ [12], [L3| that the 
uncertainties resulting from the incompleteness of this information are entirely negligible 
at x > 0.05. Moreover, these uncertainties are small even at much lower x, down to 
x ~ 10 -4 at not too small scales, Q 2 > 10 GeV 2 fll3| . Thus analyses of structure functions 
in DIS (and of total cross sections for Drell-Yan lepton-pair production JTJJ) can be 
promoted to NNLO over a wide kinematic region. Besides more accurate determinations of 
the parton densities, such analyses facilitate a considerably improved theoretical accuracy 
A t ha s (M|) ~ 0.002 of the determinations of the strong coupling. 



In the present article we extend, for x > 10~ 2 , our treatment |ll| of the flavour non- 
singlet (NS) sector dominating a^-extractions from DIS to the next-to-next-to-next-to- 
leading order (N 3 LO). This extension is facilitated by two circumstances: The first is 
the existence of constraints on the three-loop coefficient functions which prove to be 
sufficiently restrictive in this region of x. The seven lowest even- integer and odd- integer 
moments have been computed [0, |j for the structure function F 2j ns m electromagnetic 
DIS and F^ +u in charged-current DIS, respectively. Furthermore the four leading large-x 



terms of these functions are known from the soft-gluon resummation [15|, [Uj . The second 
circumstance is the rapid convergence of the splitting-function expansion in the usual MS 
factorization scheme also employed in refs. M, §]. Already the impact of the three- loop 
splitting functions is small at x>10~ 2 , in absolute size (less than 1% on a s (M|) ) as well 
as compared to the two-loop coefficient functions |ll], [12| . Hence one can safely expect 
that the effect of the unknown four-loop splitting functions on determinations of a s (M|) 
will be well below the 1% accuracy we are aiming at. 

1 See refs. || [Ij for NNLO analyses based on the integer Mellin-iV results of refs. || only. 
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As demonstrated below, the N 3 LO approximation suffices for achieving this accuracy 
in the region x < 0.75 usually covered by analyses of DIS data Terms beyond this order 
are relevant at x > 0.8, on the other hand, mainly due to the presence of large soft-gluon 
logarithms up to [ln 2Z_1 (l— x) ]/ (1—x) in the /-loop coefficient functions. The resummation 
of these logarithms (T5|, [16| has been extended to the next-to-next-to-leading logarithmic 
(NNLL) accuracy recently ||17|1 . Here we will study the predictions of this resummation for 
the factorization-scheme independent ('physical') kernels governing the scaling violations 
('evolution') of the non-singlet structure functions. Other approaches to estimate higher- 
order corrections to these kernels, not restricted to very large x, include Pade summations 



of the perturbation series [18] as well as renormalization scheme optimizations such as the 



principle of minimal sensitivity (PMS) and the effective charge (ECH) method . 



We will also compare these estimates to the full NNLO and N 3 LO evolution kernels, and 
investigate the resulting predictions at order a b s (N 4 LO) and beyond. 

The outline of this article is as follows: In section 2 we express the physical evolution 
kernels, up to N 4 LO and NNLL accuracy, in terms of the corresponding splitting functions 
and coefficient functions. The information on the three-loop coefficient functions for 
-P2,ns an d F 3 discussed above is employed in section 3 to derive approximate expressions 
for their x-dependence. Besides these functions the N 3 LO evolution kernels also involve 
convolutions of lower-order coefficient functions for which we provide compact expressions 
in section 4. These results are put together in section 5 to study the effects of the N 3 LO 
terms on the evolution of the structure functions and on the resulting determinations of 
the strong coupling constant. In section 6 we discuss the predictions of the soft-gluon 
resummation and of the Pade, PMS and ECH approximations. Finally our results are 
summarized in section 7. Some relations for the convolutions in section 4 and for the 
Pade approximations in section 6 can be found in the appendix. 
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2 Fixed-order and resummed evolution kernels 



For the choice fi 2 = fi 2 = Q 2 of the renormalization and mass-factorization scales, the 



structure functions 

F\ — 2i*i,NS j 

are in perturbative QCD given by 



^2 



X 



2,NS 



Fa 



1 v±v 



(2.1) 



F a (x,Q 2 ) = [C a {Q 2 ] ® q a ^{Q 2 )]{x) 

= E«iW 2 ) kz®g a ,Ns(Q 2 )lW 



1=0 



1=0 Jx y y y 



(2.2) 



Here c a j represents the Z-loop non-singlet coefficient functions with c a ,o(x) = 8(1 — x), and 
q a ,NS stands for the respective combinations of the quark densities. The scale dependence 
of the running coupling of QCD, in this article normalized as 



a, = 



An 



is governed by 



da q 



(5(a s 



(2.3) 



(2.4) 



dlnQ 2 l=Q 

Besides /?o and /3i [0 also the coefficients /?2 and (3% have been computed |22| in the 
MS renormalization scheme adopted throughout this study. All these four coefficients, 



A) 
Pi 

fa 



11 - 2/3 N f 
102 - 38/3 N f 

2857/2 - 5033/18 iV^ 325/54 N 2 f 
29243.0 - 6946.30 N f + 405.089 N 2 f + 1093/729 N 



f > 



(2.5) 



are required for N 3 LO calculations. The irrational coefficients in Eq. ( |2.5| ) and in Eqs. 
( p,15| ) and ( |2.18| ) below have been truncated to six digits for brevity. Nf denotes the 
number of effectively massless flavours (mass effects are not considered in this article). 
Finally the evolution equations for the parton densities in Eq. (2^2) read 

d 



dlnQ- 



[X 



Eai +1 [P^®g a , NS (Q 2 )l(x) 



1=0 



(2.6) 



where ® abbreviates the Mellin convolution written out in the third line of Eq. 
Like the coefficient functions c at i(x), the (Z + l)-loop splitting functions P a ,i(x) are scale 
independent for the above choice of /i r and \if. 
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Explicit expressions up to order a 2 can be found in refs. |23| and |J for the non-singlet 
splitting functions and coefficient functions, respectively. For the third-order splitting 



functions P a ,2(x) we will employ our approximate expressions of ref. [13]. The three-loop 
coefficient functions c a ^(x) are the subject of section 3 below. 

It is convenient to express the scaling violations of the non-singlet structure functions in 
terms of these structure functions themselves, thus explicitly eliminating any dependence 
on the factorization scheme and the scale /i/. The corresponding 'physical' evolution 
kernels K a for ft 2 = Q 2 can be derived by differentiating Eq. ( |2.2| ) with respect to Q 2 by 
means of Eqs. ( |2.4| ) and Q2.6|), and finally eliminating g aNS using the inverse of Eq. (|2.2|). 
Suppressing the dependences on x and Q 2 one arrives at the evolution equations 

Fa = [Pa + P — ;— 5 ® Ta 



dhiQ 2 \ da s 



1=0 



|a s P + E a i +1 [Pa,l-Y,PkC a .. 
I 1=1 \ fc=0 



l-k 




(2.7) 



with 



C a ,l = C a ,i 

C a ,2 — ^ C a ,2 — C a l 

C a ,3 = 3 C aj3 - 3 C a , 2 <8> C a ,l + C®J (2.8) 

C a ,4 = 4 C a , 4 - 4 C a ,3 <g> C 0)1 - 2 C®| + 4 C a , 2 <g> eg 2 - eg . 

In Eq. (|2.8|) we have used the abbreviation f® 1 for the (Z— l)-fold convolution of a function 
/(x) with itself, i.e., f® 2 = f ® / etc. The generalizations /C aj / of the kernels K a> i in 
Eq. (|2.7| ) to fi 2 7^ Q 2 can be obtained by expanding a s (Q 2 ) in terms of a s (fi 2 ,) and 
L = ln(Q 2 //i 2 ) , yielding 



fca,0 


= ^a,0 






£a,l 


= #a,l 


— (3 LK a:0 




^a,2 


= ^a,2 


- 2p LK a)1 - 


(/3iL - (3qL 2 ) K afi 


^a,3 


= ^a,3 


- 3p LK a , 2 - 


{2(3 X L - 3(3 2 L 2 ) K aA 




-o 


3 2 L - 5ftA,2/ 


1 + $L 3 ) K afi 


/Ca,4 




- AfoLK^ - 


(3p 1 L-Q(3 2 L 2 )K a , 2 



- {2(3 2 L - l^L 2 + 4/3 3 L 3 ) K a>1 (2.9) 

- (foL - 3(3 2 p L 2 - \0{L 2 + H foffiL 3 - /3 4 L 4 ) K afi . 

2 Up to the fifth order this expansion can be read off from the K a $ terms of Eq. J2.9| ). 
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At N m LO the terms up to I = m are included in Eq (|2.7|) . Hence Eqs. ( f2~8|) and (|2.9|) 
formally specify the evolution kernels K, a up to N 4 LO. Their extension to higher orders 
is straightforward but irrelevant for the time being, as at least the coefficient functions 
beyond four loops will not be calculated in the foreseeable future. 

The leading terms of the coefficient functions for x — > 1, however, are known to all 
orders from the soft-gluon resummation [15], |16, |17]. Switching to Mellin moments, 



./ 



N 



dxx N - l f{x) 



(2.10) 



for the remainder of this section, the large- N (large-x) behaviour of the coefficient func- 
tions in Eq. (|2.2|) takes the form 



CL = g (a s )e W {]nN gi (X)+g 2 (X) + a s g 3 (X)+O(a 2 f(X))} (2.11) 

up to terms which vanish for N^oo. Here we have used the abbreviation 

a* 



X = a s /3 \nN 



Air 



(2.12) 
Eq. dTT]) 

leads to the following expression for the resummed kernel up to next-to-next-to-leading 
logarithmic (NNLL) accuracy [|17 



and we have again put ji 2 r 



Q 2 . By virtue of the first line of 



N 



- (At a s + A 2 a 2 s + A 3 a 3 s ) In N 



1 + tt^s + —at A — — 



dX 



d 



- (a s p + a%)X^ - a 2 J -^(Xg 3 (X)) + 0(a 3 a (f(X)) 



dX 



(2.13) 



Thus the leading logarithmic (LL), next-to-leading logarithmic (NLL) and NNLL large- N 
contributions to the physical evolution kernels are of the form (a s In iV) n , a s (a s In N) n and 
a 2 (a s In N) n , respectively. This is in contrast to the coefficient functions which receive 
contributions up to (a. In 2 N) n . The constants A { in Eq. (|2~T3|) are the coefficients of the 

of the /-loop MS splitting functions — recall that 



x 



leading [24j] large-x terms 1/[1 

f N = -lnJV + 0(l) for f{x) 



1/[1 



x\ 



(2.14) 



As in Eq. ( p.5|) inserting the numerical values for the QCD colour factors Ca and Cp, 
these constants are given by 



Ax = 16/3, A 2 = 66.4732- 160/27 N f , 
and the yet approximate, but sufficiently accurate three-loop result [§, |3 
A 3 = (1178.8 ± 11.5) -(183.95 ±0.85)^-64/81^ . 



(2.15) 



(2.16) 
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Inserting the explicit form of the functions g l , g 2 [ 15 ] and g 3 [|T7[] in Eq. ( |2.13| ) and restoring 
the dependence on L = ln(Q 2 //U 2 ) leads to 



K 



N 



ln(l-A) + a s LA 1T ^- 

PO l — A 



a 2 B L 2 Atf Q 



A(A 



2(1 - A) 2 



+ a s [Axfc (A + ln(l - A)) + (B 1 - A lle )[3 2 \ - A 2 f3 \} 



[3 2 (1 - A) 



+ a 2 s L 



[B L - 7e Ai)#, - A 2 



A(A-2) 



Alp! 



ln(l - A) 
- A) 2 



^i(Te + Ca)0o + 2A 27e /5 + A 3 - 2E l7e/ 3 2 - 25 2/ 3 



A(A-2) 
2/5 (l-A) 2 



+ 



+ 2D 



2(A l7e - 50A/5 2 ln(l - A) + A x 0{ (A 2 - ln 2 (l - A) 

• ^i&A) A 2 - Az/hPo (A 2 - 2A - 2 ln(l - A) 
A(l-A) 



2/? 3 (l " A) 1 



[l-2\y 



+ 0(af/(A,L)) 



(2.17) 



— ^res,Z 



Z=0 



Here (2 = k 2 /6, and 7e represents the Euler-Mascheroni constant, 7e ~ 0.577216. Fur- 
thermore Bi = —4 [|15j, and the constants B 2 and D 2 are related by fIT| 



B 2 + D 2 = 36.2657 + 6.34888 A7 . 
We will return to the latter coefficients at the end of the next section. 

m+l 



(2.18) 



After subtracting the terms up to order a™ in /C^ s already taken into account in 



the N m LO terms (|2.9|) , Eq. (|2.17|) can be added to these fixed-order results to obtain the 
(N m LO + resummed) approximation for the non-singlet evolution kernels, 



K~ N — \ n l + l fk" N — k" N \ _1_ \c 

^N^LO+res — a s K^a,! ^resjJ "r /V 



iV 
res 



(2.19) 



1=0 



Due to the renormalon singularities at A = 1 and A = 1/2 in Eq. (|2.17 ) the resummed 
evolution equations cannot be uniquely inverted to x-space, unlike the fixed-order approx- 
imations tCaiJ-Q. Note that strength of these singularities — located at iV ~ 2000 and 
iV ~ 45 for A = 1 and A = 1/2, respectively, at a s = 0.2 and Nf = 4 — increases with the 
order of the soft-gluon expansion: the behaviour is logarithmic at the leading-log level, 
but involves poles of order k for the N fc LL approximations. For our numerical study of 
the all-order case at the end of section 6 we will use the standard 'minimal prescription' 
contour [27] for the Mellin inversion. This contour runs to the left of the renormalon 
singularities, but to the right of all other poles in the iV-plane. 
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3 The 3-loop non-singlet coefficient functions 



The /-loop coefficient functions c a j for the non-singlet structure functions J- a =i,2,3 defined 
in Eq. ( |2.1| ) can be represented as 

21-1 

c a M N f) = E A M V m + A S(l -x)+ th (x, Nf) 

m=0 
21-1 

+ £K?^ + ^o)- (3.1) 

n=l 

Here A\ , B h C^j and D^j are numerical coefficients which in general depend on the 
number of flavours Nf, and we have employed the abbreviations 

"ln fc (l - x) 



x 



L x = ln(l -x) , L = In 2 (3.2) 



for the -(—distributions (see Eqs. ( |3.3|) and (|3]4]) below) and the end-point logarithms. 
The functions c^f ooth (x, Nf) collect all contributions which are finite for < 2 < 1. This 
regular term constitutes the mathematically complicated part of Eq. (|3.1|) , it involves 
higher transcendental functions like the harmonic polylogarithms introduced in ref. |28| . 
As usual, the -(--distributions are defined via 

dxa(x) + f(x) = / dx a{x) {f {x) — f {!)} (3.3) 
Jo 

where f(x) is a regular function. The convolutions with the distributions occurring in 
Eq. ( |3.1p can be written asQ 

+ xf(x)-r^—\n t+ \l-x) . (3.4) 
k + 1 

As already indicated in Eq. (|3~T|), the coefficients of V m and of 5(1 — x) are independent 
of the choice of the structure function. 

The three-loop contributions c s 3 known from the soft-gluon resummation read 

c S3 (x,3) = 512/27 (V 5 -Ll) -14400/81 P 4 + 264.062 V 3 + 1781.704 V 2 
c 5>3 (x,4) = 512/27 (£> 5 - ^i) - 13760 / 81 £> 4 + 188. 210 £> 3 + 1962. 178 £> 2 (3.5) 
c s>3 (x,5) = 512/27 (V b -L\) -13120/81 £> 4 + 113.938 £> 3 + 2131.195 V 2 , 

where we have again truncated the irrational coefficients and restricted ourselves to the 
practically relevant cases Nf = 3, 4 and 5. Besides the P m -terms determined in ref. fT6|| . 

3 The second line of Eq. ( |3.4| ) is given by —xf(x) J^dya(y) for a general +-distribution [a(x)]_|_. 
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Eq. ( |3.5| ) also includes the leading integrable large-re logarithm. The general relation 
between the coefficients of this term and the leading -(—distribution has been conjectured 
in ref. [p9 |. Eqs. ( |3.5|) complement the main present constraints on c a 3 {x, Nf) provided by 



the computation 0, || of the seven lowest even- integer and odd-integer moments ( [2.10|) , 
respectively, for electromagnetic (e.m.) DIS and the charged-current (CC) combination 
F 3 U+U . Note that the coefficients of the leading small- a; logarithms are presently unknown 
here, unlike for the splitting functions and the singlet coefficient functions [TIJ . 



We use this information for approximate reconstructions of c 2 3 (x, Nf) and c 3 3 (x, Nf) 
at Nf = 3, 4 and 5. Our method is analogous to the treatment of the three-loop splitting 
functions in refs. [[H], [12], 13]: A simple ansatz is chosen for c^ ooth in Eq. ( |3.1|) , and its 



free parameters are determined from the available moments together with a reasonably 
balanced subset of the coefficients A 3 °' 1 ' ) , B 3 , cj^ and D^ 3 \ The ansatz for c^™ ooth and 
the choice of the non- vanishing end-point parameters are then varied in order to estimate 
the residual uncertainty of c a3 . Specifically we keep A 3 ; one of each pair 

4 0) and B 3 , 

C^l and C^ 3 3 , C^l and C^\\ one or two of the D^ 4 ^; and one or two parameters of 
a polynomial up to second order in x representing c^ ooth . For a few combinations the 
resulting system of linear equations which fixes these parameters by the seven moments 
becomes almost singular, resulting in exceptionally large numerical coefficients. After 
rejecting those about 5% of the combinations for which the modulus of at least one 
parameter exceeds 10 5 , we are left with about 90 approximations for each case. 

Before we present the approximate results for c a3 (x, Nf), it is appropriate to illustrate 
our procedure by applying it to a known quantity, for which we choose the two-loop e.m. 
coefficient function c 22 (x,Nf = 4). Adopting the coefficients of T> 3 and T> 2 defined in 
Eq. ( |3.2|) from the soft-gluon exponentiation, the procedure described in the preceding 
paragraph is applied to this function with the small adjustment that two of the C 2 2 ' 
are kept as does not occur at two loops according to Eq. ( |3.1|) . Also here we reject a 
couple of combinations, those with parameter(s) of modulus above 3 • 10 3 . The remaining 
about 70 approximations are compared to the exact result of ref. [f| in Fig. 1. 

The seven lowest even-integer moments supplemented by the soft-gluon coefficients 
^4 2 m>1 - ) prove to constrain c 22 (x) rather tightly at x>0.3. The region 0.1<x<0.3 is 
less accurately covered, and at x< 0.1 the lack of small- x information mentioned above 
becomes very prominent. Also shown in Fig. 1 is the exact result || for c 22 in charged- 
current DIS. The difference to the electromagnetic case — originating in a sign difference 
of the contributions from 'y/W +q^q+q+q with identical quarks in the final-state 
- is clearly visible only at x < 0.2. The effect of this difference on the evolution of the 
structure functions at NNLO is unnoticeable at x > 0.1, and amounts to less that 1% for 



x > 0.01, see Fig. 11 of ref. Jllj]. We expect that the corresponding three-loop effect will 
at least not be larger. Hence the approximations for c 2 2 (x, Nf), constructed for the e.m. 
case, should be applicable also for neutrino DIS without introducing any relevant error. 
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Figure 1: Approximations for the two-loop coefficient functions c 2 2 (x, A/ = 4) for -F^ns 
obtained from the lowest seven even-integer moments and the two leading soft-gluon 
terms, compared to the exact result of ref. ||. Also shown is the corresponding exact 
coefficient function for -F 2 ns i n charged-current DIS. 



The corresponding approximations for c 2 3 and c 3 3 are shown in Fig. 2 for Nf = 4 
(concerning the scale of the ordinate recall the rather small expansion parameter ( |2.3| ) ). 
As expected, the accuracy pattern is qualitatively similar to the two-loop case of Fig. 2. 
The uncertainty of c 3 3 is smaller than that of c 2 3 at small x, since for c 3 3 the lowest 
calculated moment, N — 1, is closer to the location of the rightmost pole at N = 0. 
For both functions two representatives, denoted by A and B, are selected which rather 
completely cover the uncertainty bands. With c 53 of Eq. ( |3.5|) these representatives read 

c ^ 3 ( x ,4) = c s>3 {x, 4) - 6456.231 V x - 1085.97 5(1 -x) +258.876 L\ 

- 22430.79 Li - 74705.15 x 2 - 4062.14 (2 - x) - 313.0 L 3 

cg 3 (z,4) = c s>3 (x, 4) - 5081.227 V x + 5028.23 V + 1059.423 L\ (3.6) 

- 7292.84 L\ - 17741.28 (2 - x) - 18154.5 L + 1168.02 L 3 Q 

and 

c ^ 3 ( x? 4) = c 53 (x, 4) — 6940.648 V\ — 1526.23 V — 42.598 L\ 

- 33562.64 L x - 91639.87 x 2 - 5898.84 + 424.49 L 2 

c^ 3 (x,4) = c 53 (x, 4) — 4907.988 T>\ + 5587.906 V + 1092.436 L\ (3.7) 

- 8267.99 L 2 - 18120.78 - 7083.63 L + 283.59 Ll . 
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Figure 2: Approximations for the three-loop coefficient functions for -F 2 6 ns (left) and the 
charged- current F£ +u (right) derived from the respective seven lowest moments J7|, || and 
the soft-gluon terms (|3.5|). The full lines show the selected functions (|3.6|) and (pT 
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Figure 3: The convolution of the approximations ( |3.6| ) and ( |3.7| ) selected from the previous 
figure with a shape typical of hadronic non-singlet distributions. 
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The uncertainty bands of Fig. 2 do not directly indicate the range of applicability 
of our approximations, as the coefficient function enter the structure functions and their 
evolution only via the smoothening convolution (|2.2|) with non-perturbative initial distri- 
butions. In Fig. 3 we therefore present the convolutions of the results (|3.6| ) and ( |3.7|) with 



a typical non-singlet shape. This illustration shows that the residual uncertainties of c a 3 
do not lead to any relevant effects for x > 0.1. The situation at smaller x depends on the 
numerical size of the c a3 contributions to the evolution kernels given by Eqs. Q2.7D and 
( |2.8D . Anticipating our findings in section 5, we note that these contributions are actually 
unproblematically small for x > 10~ 2 . 

The results for Nf = 3 and Nf = 5 are similar to those presented in Fig. 2 and Fig. 3. 
For brevity they are not displayed graphically in this article. The selected approximations 
for Nf = 3 are given by 

c^ 3 (x,3) = c S3 ( X) 3) - 6973.782 V x - 3929.78 8{1 -x) + 232.676 L\ 

- 35949.75 L x - 93141. 53x 2 - 10283.51 - 418.43 
c g 3 ( X; 3) = Cs 3 ( Xi 3) - 5575.903 V x + 5474.48 V + 927.478 L\ (3.8) 

3 
J 



- 4646.38 L\ - 23345.85 - 11094.92 L + 759.69 L 3 



and 



c^ 3 (x,3) = c S3 (x, 3) -7018.496 V x - 2957.33 5(1 -x) + 162.667 L{ 

- 28363.91 L x - 74640.87a; 2 - 5720.48 (1 +x) + 330.21 L 2 

c^ 3 (x,3) = c S 3 (x, 3) — 5160.335 V\ + 6896.360 V + 1109.512 L\ (3.9) 

- 7715.43 L\ - 20541.22 - 7595.83 L + 290.34 Lq . 

The corresponding functions for Nf = 5 read 

c^ 3 (x,5) = c S3 (x, 5) - 5951.174 V x - 391.37 V - 2341.422 L\ 

+ 19986.58 L x + 5517.39 (2 - x) + 5969.63 L - 284.23 L 3 

,3 (x,5) = c s 3 (x, 5) -4802.695 V x + 3784.97 V + 1041.041 L\ (3.10) 



'2,3 V 



- 8021.15 L\ - 15556.5 (2 - x) - 16445.21 L + 1084.36 L 3 



and 



c^ 3 {x,5) = c S 3 (x, 5) — 6560.902 V x — 2412.54 V + 98.499 L\ 

- 27899.69 Lx - 82015.71 x 2 - 2983.82 - 61.43 L 3 

cj 3 (x,5) = c S3 (x, 5) -4637.854 Pi + 4317.29 £> + 1070.036 L\ (3.11) 

- 8767.685 L\ - 15676.57 - 6543.88 L + 276.32 L 3 . 



In all cases the average 1/2 (cf 3 + cf 3 ) represents our central result. 
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We conclude this section by returning to the coefficients B 2 and D2 entering the NNLL 
soft-gluon resummation of the quark coefficient functions ( [2.11| ) and structure-function 
evolution kernels (|2.17| ). If only one of these constants were present, say D 2 , then this 
constant would be fixed by the consistency of Eq. (|2.11|) with the soft-gluon part c S2 
of the NNLO coefficient functions of ref. 0], more precisely by the coefficient A 2 in 
Eq. Digressing for a moment, we note that this situation is actually realized for 

the (very closely related) NNLL soft-gluon resummations of the quark-antiquark annihi- 
lation contribution to the Drell-Yan cross section |ljj and of the cross section for Higgs 
production via gluon-gluon fusion in the heavy top-quark limit. For these two processes 
the NNLO soft- and virtual-gluon contribution have been computed in refs. |25| and f26fl , 
respectively. In the present DIS case, however, this consistency conditions only implies 
the constraint Q2.18| ). An exact result for the coefficient A 2 of the three- loop coefficient 
functions ( |3.1|) would suffice to determine B 2 and D 2 , as this coefficient is related to the 
combination B 2 + 2D 2 independent from Eq. ( |2.18| ). 



As discussed in the paragraph below that of Eq. ( |3.5|) , all our approximations (shown 
for Nf = 4 in Fig. 2) include A 2 . For about 95% of these approximations this coefficient 
falls into the range 



A 



(3) 



-6800 
-6350 
-5950 



- 5800 

- 5500 

- 5200 



for N 



f 



(3.12) 



The comparison of these results to the expansion of Eq. ( |2.11|) (using g 3 (A) of ref. |T 
leads to the rather weak constraints 



32 . 


. 87 


42 . 


. 93 


49 . 


. 98 



for N 



f 



which can to sufficient accuracy be combined to the estimate 

1 




(3.13) 



B 2 ~ -P» + - tf3 Q P« with 8 < i< 12 . 



(3.14) 



Here P l s _ 1 represent the coefficients of 5(1 — x) in the /-loop quark splitting functions. 
Retaining the colour factors Ca — N c — 3, Cf = — 1)/(2N C ) these coefficients read 



P — 3 Cp 



Pi = C|(- + 24C3-12C 2 )+C f .C^(--12C 3 + yC2) (3.15) 
for our normalization (|2.3|) of the expansion parameter. 
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4 Convolutions for the N 3 LO evolution kernels 



Besides the (Z + l)-loop splitting functions and the /-loop coefficient functions, the N'LO 
evolution kernels ( |2.7| ) involve simple and multiple convolutions of the coefficient functions 
of lower order. Required at N 3 LO are the simple and double convolutions of the one-loop 
coefficient functions c ol with themselves, cf\ and C® 3 , and the convolutions of the one- 
with the two-loop coefficient functions, c ol ®c o2 , see Eq. (|2.8|) . Especially the latter 
lead to rather complex exact expressions. These terms, however, do not require any 
attention if the evolution is carried out using the moment-space technique PD |, as in 



iV-space the convolutions reduce to products. On the other hand, many analyses of data 
on structure functions are performed using 'brute-force' x-space programs for solving the 
evolution equations. For application in such programs we provide compact and accurate 
parametrizations of the convolution contributions to the evolution kernels up to N 3 LO. 

These approximations are derived analogously to those of the two-loop coefficient 



functions in ref. The -(—distribution parts are treated exactly (truncating irrational 
coefficients), see the appendix. The integrable x< 1 terms are fitted to the exact results for 
x > 10~ 6 . Finally the coefficients of 5(1 — x) are slightly adjusted from their exact values 
using the lowest integer moments. The resulting parametrizations deviate from the exact 
results by no more than a few permille. This accuracy applies directly to Eqs. (|4.1|) 
as well as to their convolutions with typical hadronic input distributions. 



Using the abbreviations ( |3.2| ) the simple convolutions of the one-loop coefficient func- 
tions for F 2) ns and F 3 can be written as 

cfl(x) = 256/9 X? 3 -64£>2 - 283.157 £>i + 304.751 V + 346.213 5(1 - x) 



- 26.51 L\ + 192.9 L\ + 198.2 L x + 113.0 L\L, 



o 



- 1.230 Lq + 9.466 L Q + 32.45 L - 483.3 x - 410.5 (4.1) 



and 



c®?(ar) = 256/9P 3 -64P 2 - 283.157Pi + 304.751Po + 345.993 5(l-x) 

- 27.09 L\ + 162.1 L\ + 248.0 L x + 91.79 L\L 

— 1.198 Lq + 3.054 L\ + 65.54 L - 305.3 x - 335.7 . (4.2) 

The corresponding parametrizations for the double convolutions read 

cg(x) = 1024/9 X> 5 - 1280/3 V 4 - 2757.883 P 3 + 9900.585 P 2 

+ 3917.516 V x - 12573.13 V - 2851.0 5(1 - x) (4.3) 

- 151.4 L\ + 118.9 L\ - 6155 L\ - 47990 L{ - 30080 L x + 6423 L\L Q 

- 0.35 Ll - 4.30 Lg - 106.7 L 3 - 1257 L% - 4345 L + 3618 x + 8547 
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and 



cf?0) = 1024/9 V 5 - 1280/3 X> 4 - 2757.883 P 3 + 9900.585 L> 2 

+ 3917.516 V l - 12573.13 V - 2888.1 5(1 - x) (4.4) 

- 138.4 L\ + 409.0 L\ - 1479 L? - 24700 L\ + 9646 L x - 10080 L^Lq 

— 0.119 + 3.126 L^ + 84.84 L 3 + 288.7 + 264.87 L + 15410 x + 15890 . 

The convolutions c a 2 ® c a i f° r ^2,ns i n electromagnetic DIS and for the charged-current 
combination F^ +v are parametrized as 

[c 2j2 ® c 2jl ](x) = 

1536/27 P 5 - 343.702 P 4 - 633.29 P 3 + 5958.86 P 2 - 6805.10 ^ - 2464.47 V 

- 101.7 L^ - 155.1 L\ - 6553 - 23590 L\ - 10620 L x + 9290 L^L - 0.35 L 5 
+ 0.64 Lq + 92.93 Lq + 761.9 Lg + 2450 L - 1251 x + 6286 + 8609.2 5(1 - x) 

+ N f {7.0912 L> 4 - 55.3087X? 3 + 18.629 V 2 + 619.865 V x - 584.260 V 

- 11.71 + 60.82 - 618.0 L^ - 1979 L x - 919.6 L\L Q + 0.48 L 4 (4.5) 

- 1.08 Lq 1 - 43.83 Ll - 125.5 L - 295.1 x + 522.4 - 809.145(1 - a:)} 

and 

[ C 3,2 ® C 3)1 ](x) = 

1536/27 £> 5 - 343.702 £> 4 - 633.29 £> 3 + 5958.86 P 2 - 6805.101?! - 2464.47X? 

- 77.39 L\ + 289.1 L\ - 2823 Lf - 12500 L\ + 25420 Lx + 9515 L\L - 0.524 L 5 

- 6.104 Lq + 39.23 Lq + 553.5 Lq + 1393 L + 20080 x + 2548 + 8478.2 5(1 - a;) 

+ N f {7.0912^4 - 55.3087 L> 3 + 18.629 L> 2 + 619.8651?! - 584.260X? 

- 14.30 L\ + 10.47 L\ - 775.2 L\ - 2458 L x - 392.9 L\L Q + 0.482 L 4 (4.6) 
+ 2.541 L\ - 41.04.Ljj- 223.9 L - 891.9 x + 468.0 - 803.43 5(1 - x)} . 

The corresponding results for the charged-current quantities i^NS an d ^a" - " are obtained 
by replacing the second and third line of Eqs. (|4.5| ) and ( 14.61) , respectively, by 

- 109.2 L\ - 243.4 L\ - 6890 L\ - 24000 L\ + 10840 L x + 9144 LjL - 0.45 L 5 
+ 1.80 Lq + 114.0 L 3 + 856.6 Lq + 2602 L - 711.6 x + 6298 + 8569.2 5(1 - x) 

and 

- 77.39 L\ + 295.5 L^ - 2587 L\ - 10580 Lj + 30580 L x + 6461 L^L - 0.404 Lq 

- 5.525 Lq + 23.80 L 3 + 484.7 Lq + 1577 L + 22220 x + 3349 + 8485.0 5(1 - x) . 
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5 Numerical results for the scaling violations 



In this section we illustrate the effect of the next-to-next-to-next-to-leading order (N 3 LO) 
contributions to the physical evolution kernels ( |2 . 7| ) — (pT9|) for the electromagnetic structure 
function -F 2 ns an d the charged-current combination F^ henceforth simply denoted F. 



3- 



Specifically, we will discuss the logarithmic derivatives F a = d\iaF a /d\iaQ 2 calculated at 
a fixed reference scale Q 2 = Ql for the initial conditions 

F 2>m ( x >Ql) = ^zMl) = x°-\l-x) 3 . (5.1) 



The simple model shape ( |5.1|) incorporates the most important features of non-singlet 
x-distributions of nucleons. Its overall normalization is irrelevant for the logarithmic 
derivatives considered here. The reference scale Ql is specified via 

a s (tf = Q 2 ) = 0.2 (5.2) 

irrespective of the order of the expansion. Eq. ( |5.2j ) corresponds to Ql ~ 30 GeV 2 , a scale 
typical for fixed-target DIS, for a s (M|) ~ 0.116 beyond leading order (LO). The same 



input (|5.1|) and ( |5.2| ) is chosen at all orders and for both structure functions in order to 
facilitate a direct comparison of the various contributions to the evolution kernels ( [2.7|) . 
All graphical illustrations below refer to Nf = 4 effectively massless quark flavours. 

Before turning to the numerical results we have to specify our treatment of the four- 
loop splitting functions P a 3 . These functions enter Eq. ( |2~7| ) at order a\ (N 3 LO) together 
with the three-loop coefficient functions ( |3.6| ) - ( |3.11| ) and the convolutions ( ^4.3| ) - ( |4.6|) . 
As already mentioned in the introduction, the size of the two- and three-loop terms in 
the expansion of the non-singlet splitting functions strongly indicates that the effects of 
P Qj 3 are very small in the x-region addressed by the present study, x > 10 -2 . Hence a 
rather rough estimate of these quantities is sufficient here. We have checked that the 
[0/1] Pade approximation^] gives a reasonable, though not particularly accurate estimate 
of the three- loop non-singlet splitting functions P^ 2 i n A^-space for Nf = 3 ... 5. Thus we 
choose the Mellin inverse of 

P a N s - V [P%1 [i/i] Fade , V = 0...2 (5.3) 

as our estimate of P a>3 (x), i.e., we assign a 100% error to the predictions of the [1/1] 
Pade summation (the [0/2] Pade results are similar). The results obtained by combining 
Eq. (|5.3| ) for r\ = 2 with c^ 3 in Eqs. (|3.6|) - ( |3.11| ) are denoted by N 3 LO J 4 in the figures 
below, those using 77 = and cf 3 by N 3 LOb- As in section 3 the central predictions 
1/2 (N^LO^ + N 3 LO B ) are not shown separately. 



The logarithmic scale derivatives F 2 ns an d F 3 resulting from Eqs. (|5\1 



are 
2 



shown in the left parts of Fig. 4 and Fig. 5, respectively, for the standard choice fif = Q 

4 A brief discussion of the Pade summations and the resulting higher-order approximations can be 
found in the next section. 



15 




0.03 



0.02 



0.01 



-0.01 



-0.02 



-0.03 



Figure 4: The perturbative expansion of the scale derivative -F 2 NS = d\nF 2 ^s/dlia.Q 2 
of the electromagnetic structure function F 2j ns at fx 2 = Q 2 ~ 30 GeV 2 for the initial 
conditions specified in Eqs. (|5.1| ) and ( |5.2|) . The differences between the predictions at 
different orders in a s are shown on a larger scale in the right part. 
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Figure 5: As Fig. 4, but for the charged-current combination F s = F^ +v . In all figures 
the subscripts A and B at N 3 LO refer to the approximations discussed below Eq. 
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of the renormalization scale. At x £ 0.5 the size of the NNLO and N 3 LO corrections 
can barely be read on the scale of these graphs, therefore the differences F a — (F a ) nnlo 
are displayed on a larger scale in the left parts of both figures. The difference of the 
N 3 LO A and N 3 LO B results is very small down to x ~ 10 even on this enlarged scale, 
demonstrating that our approximations (3J3) - ( 3.11 ) and (|5.3|) are completely sufficient 



in this region of x. For both structure functions the N 3 LO corrections are large only at 
very large x, where the kernels are dominated by the universal soft-gluon contributions. 
Towards smaller x the N 3 LO effects rapidly decrease, e.g., from 6% at x — 0.85 to 2% 
at x — 0.65. The corresponding NNLO contributions amount to about 12% and 6%, 
respectively. At 10~ 2 < x < 0.6 the N 3 LO corrections are particularly small for -F 2j ns- 
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Figure 6: The dependence of the NLO, NNLO and N 3 LO predictions for cHni^Ns/^hiQ 2 
at Q 2 = Ql ~ 30 GeV 2 on the renormalization scale [i r for six typical values of x. 
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Figure 7: The renormalization scale uncertainties, as estimated by the quantities Aab s F a 
defined in Eq. fl5.4Q , of the perturbative predictions for the scale derivatives of F 2t NS 
and F 3 displayed in Fig. 4 and Fig. 5, respectively. Also included are the corresponding 
approximate N 4 LO results derived in section 6 below. 



The dependence of these scale derivatives on the renormalization scaleQ is illustrated 
in Fig. 6 and Fig. 7. In the former figure the consequences of varying \i r are shown for 
-p2,NS at six representative values of x (note that the scales of the ordinates are different in 
all six parts). Here we vary fi r over a rather wide range, | Q 2 < fi 2 < 8Q 2 , corresponding 
to 0.29 > a s (Hr) ~ 0-15 for the initial condition ( |5.2|) . In the latter figure we display the 
absolute scale uncertainties of F 2 ns an d F 3 at Q 2 = Ql, estimated by 



'0' 

. _ 1 r 1 _ 



l - {max [PM = \q 2 -.. 4Q 2 )] - min [P a {£ = 1 Q 2 . . . 4Q 2 )] } , (5.4) 

i.e., using the smaller conventional interval \ Q . . .1Q for fi r . Also shown here are the 
further improvements resulting from including the approximate a 5 s (N 4 LO) contributions 
to Eqs. ( p.7|) - ( ^9|) discussed in the next section. 



Our new N 3 LO results represent a clear improvement over the NNLO stability JTT| for 
all x-values of Fig. 6 except for x = 0.05 (here, however, the absolute spread is very small, 
see Fig. 7), where the difference between the N 3 LOa and N 3 LOb results at small ji r be- 
comes comparable to the /i r - variation at NNLO. This enhanced sensitivity at small scales 
is due to the larger values of ot s up to almost 0.3, which enter the approximate contribu- 
tions to Eq. (|2.7| ) as a 4 . The present approximation uncertainties of the N 3 LO results are 



As already mentioned below Eq. (|2.4|), we use the MS renormalization scheme throughout this study. 
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actually dominated by the estimate ( |5.3| ) for the four-loop splitting functions, not by the 
residual uncertainties of the three-loop coefficient functions quantified in Eqs. (fT6| ) - (|3.1i|) . 

The results shown in Fig. 7 correspond to relative uncertainties (AabsF a )/F a of 8% 
at NNLO and 5% and 3% at N 3 LO and N 4 LO, respectively, for both F 2iNS and F 3 at 
x = 0.85. The corresponding figures at x = 0.65 read 5% (NNLO), 2% (N 3 LO) and 1% 
(N 4 LO). These scale uncertainties are rather similar to the relative size of the highest- 
order contributions at /z 2 . = Q 2 in Fig. 4 and Fig. 5 (for the N 4 LO contribution see Fig. 11 
in section 6). Hence the \x r - variation ( |5.4| ) and the size of the last contribution included 
in Eq. Q2.7| ) yield consistent uncertainty estimates in this region of x. At smaller x the 
absolute scale uncertainties are very small at N 3 LO and N 4 LO. For A abs F 2N s values even 
below 0.001 are reached for x < 0.5 at N 3 LO and x < 0.6 at N 4 LO. 

We conclude this section by illustrating the effect of the higher-order terms in Eq. fl2.7|) 
on the determination of a s from the scaling violations of non-singlet structure functions. 
For this illustration we assume that F 2j ns and F 3 at Ql ~ 30 GeV 2 are given by Eq. ( |5.1| ) 
with negligible uncertainty. The resulting average N 3 LO predictions of i^.NS and F 3 for 
fi 2 = Qq and a s = 0.2 are employed as model data at Xk = 0.1 k — 0.05 with k = 1, . . . ,8. 
Roughly following the experimental pattern, we assign errors of 0.005 for k — 2, ... 6, of 
0.01 for k = 1, 7 and of 0.02 for k = 8 to these data points (for Eq. (|5.5| ) only the relative 
size of these errors is relevant). Again already including the N 4 LO estimate obtained in 
the next section, the fits of a s (Ql) to these model data yield 





)nlo = 


0.2080 


+ 0.021 

- 0.013 ' 


0.2035 


+ 0.019 
- 0.011 




'nnlo = 


0.2010 


+ 0.008 
- 0.0025 ' 


0.1995 


+ 0.0065 
- 0.0015 




'n 3 lo = 


0.2000 


+ 0.003 
- 0.001 ' 


0.2000 


+ 0.0025 
- 0.0005 




'n 4 lo = 


0.2000 


+ 0.0015 
- 0.0005 ' 


0.2005 


+ 0.0015 
- 0.0005 



where the first column refers to F 2 ^s an d the second to F 3 . The central values represent 
the respective results for fi 2 = Ql, and the errors are due to the renormalization scale 
variation | Q 2 < fi 2 , < 4Q 2 , for the N 3 LO and N 4 LO cases combined with the approxima- 
tion uncertainties. Unlike the NNLO terms, the N 3 LO and N 4 LO corrections do not cause 
significant shifts of the central values, but just lead to a reduction of the /i r uncertainties 
which reach about ±1% at N 3 LO. The difference of the NLO and NNLO central results 
for F 3 is half as large as that for i^.NS- This effect is due to larger positive corrections 
to the logarithmic derivative at x < 0.4 in the former case (see Figs. 4 and 5), which 
counteract the effect of the negative corrections at large x in the fit. As far as Eq. (|OD 
can be compared to the fits of real data in refs. 0, [| (where higher-twist contributions 
affecting the central values are included), our finding are consistent with those results. 
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6 Resummations and optimizations 



Finally we address the predictions of the soft-gluon resummation, the ECH and PMS 
scheme optimizations, and the Pade approximations for the physical evolution kernel K a 
in Eq. (|2.7|) . Being especially interested in the region of large- x / large- N, where the 
higher-order corrections are large but similar for -F 2j ns an d F 3 , we will for brevity focus on 
the former, more accurately measured quantity. Also in this section the numerical results 
are given for Nf = 4 and the initial conditions ( |5.1| ) and ( |5.2| ). We will mainly consider 
the predictions of the above-mentioned approaches at fixed order in a s , and only at the 
end briefly turn to the all-order results for the soft-gluon exponentiation and the Pade 
summations. 



The N LO predictions of the soft-gluon resummation for the kernels (|2.7|) are given by 
the terms a z s +1 /C res ^ in Eq. ( |2.17|) . Recall that the leading logarithmic (LL), next-to- leading 



logarithmic (NLL) and next-to-next-to-leading logarithmic (NNLL) contributions behave 
as ln' +1 iV, ln'iV and (at l> 2) ln^iV, respectively. The terms in the /-loop coefficient 
functions c a ,i proportional to ln fc iV with k = I + 2, ... ,21 cancel in the combinations 
( |2.8| ) for I > 2. This implies that, from I = 5 onwards, actually none of the four leading 
ln fc iV terms of c a ,i presently fixed by the soft-gluon exponentiation (|2.11|) contributes to 
the N'LO kernels (|2.7| ). Consequently, we expect a pattern for the numerical soft-gluon 
approximations to the physical kernels which is rather different from that discussed in 



ref. for the MS coefficient functions. 



The cumulative effect of soft-gluon terms at NNLO (a 3 ) and N 3 LO (a 4 ) is compared 
to the (approximate) full results in Fig. 8. The two solid curves in the right plot refer to 
the N 3 LOa and N 3 LOs approximation discussed below Eq. (|5.3| ), the two NNLL results 
to £ = 8 and £ = 12 in Eq. (|3.14j ). For the (undisplayed) NLO contribution the LL and 



NNLL predictions are considerably smaller and larger, respectively, than the full result, 
whereas the inclusion of also the a^N term arising from soft and virtual emissions leads to 
a reasonable approximation. Combined with this situation, the results of Fig. 8 indicate 
that the number of soft-gluon logarithms required for a realistic approximation at N'LO 
systematically increases with the order /: The full NLO, NNLO, and N 3 LO curves run 
between the LL and NLL, close to the NLL, and between the NLL and NNLL results, 
respectively. The NNLL soft-gluon contribution may thus be expected to represent a 
reasonable estimate for the N 4 LO (a^) term of Eq. ( |2.7| ) at large x / large N. As shown 
in Fig. 90, however, the spread due to the present uncertainty ( p. 14 ) of the parameters 



B2 and D2 entering Eq. (|2.17|) is unfortunately rather large. Moreover, even if with 
this uncertainty removed, e.g., by a future exact calculation of the three- loop coefficient 
functions, Figs. 8 and 9 indicate that the soft-gluon resummation can hardly be expected 
to provide accurate information on the N 4 LO term, even for moments as large as N ~ 30. 

6 Already Fig. 3 demonstrates that the universal soft-gluon terms do not provide a good approximation 
for x < 0.7, hence the corresponding a;-space results in Figs. 9 and 12 are shown only at larger x. 
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Figure 8: The successive soft-gluon approximations for the NNLO (left) and N 3 LO (right) 
contributions to the moments (|2.10| ) of the evolution kernel ( ^7|) at fi 2 , = Q 2 , compared 
with the full results for i^NS addressed in the previous section. Besides the ln n N terms 
of K reS 2 in Eq. (2.17), also the N° contribution is included for the NNLO case. 
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Figure 9: The results of the soft-gluon resummation (|2.17| ) for the N 4 LO contribution to 
the kernel ( |2.7| ). The left part corresponds to the right plot of Fig. 8, in the right part 
the resulting large- a; predictions are shown for the corrections to the results of Fig. 4. 
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The renormalization-scheme optimizations assume that that the higher-order correc- 
tions to the N'LO physical kernels Kjp N in iV-space given by 

= K®* = a s K% (1 + a s rl x + ... + (6.1) 

are small in a certain 'optimal' scheme. The principle of minimal sensitivity (PMS) 
proposed in ref. |19|] selects this scheme by the requirement 



~dWT = ' ( } 

where d/d(RS) abbreviates the derivatives with respect to the I independent parameters 
specifying the renormalization scheme at N'LO. In the effective charge (ECH) method of 
ref. [p0| , on the other hand, these parameters are chosen such that 



<i = ••• = <i = • (6.3) 

Assuming that in these schemes the next terms r^ l+l are not just small but vanishing, 
the transformation back to MS (or any other scheme) leads to the respective PMS and 
ECH predictions for this quantity in terms of r~ \ . . . r^i and the coefficients ( j2.4| ) of the 
/3-function. Up to r 4 these predictions are explictly given in Eqs. (6) - (11) and (13) - (17) 
of ref. f31| , thus we refrain from repeating them here. 



Another approach for estimating the higher-order corrections is provided by the Pade 



summation of the perturbation series, for QCD in detail discussed, e.g., in refs. ||18|| . In 
this method K.%' N in Eq. ( |6.1| ) is replaced by 



^a,[N/V\ ~ a « A 1 N v N (0.4) 



with 

V > 1 and Af + V = I . (6.5) 

The determination of the parameters pt and qj from the r% . . . ri of Eq. (|6.1|) are autom- 
atized in programs for symbolic manipulation such as MAPLE |32] . Expanding K. 1 ^ uf/D] 
to order Z + l then yields the [N"/T>] Pade predictions for the N l+1 LO coefficients r^ l+1 . 
Also these predictions need not to be written down here. Beyond the second-order results 
there is no obvious relation between the predictions of the scheme optimizations and those 
of the Pade approximations. Consistent result of these methods for r; >2 are thus usually 



considered as evidence of the approximate correctness of these predictions [18 



The PMS, ECH and Pade results for the NNLO and N 3 LO iV-space kernels (0) are 
compared in Fig. 10 to the (approximate) full results already shown, on the same scale, in 
Fig. 8. Disregarding large relative, but small absolute deviations at NNLO for small N, 
the PMS and ECH results (which are very similar at NNLO and identical at N 3 LO) 
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Figure 10: The PMS, ECH and Pade estimates of the NNLO (left) and N 3 LO (right) 
parts of the iV-space evolution kernel for F 2j ns at y? r = Q 2 , compared with the full results 
illustrated in x-space in section 5. The scales of the graphs are the same as in Fig. 8. 
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Figure 11: The PMS, ECH and Pade estimates of the N 4 LO contributions to the evolution 
kernel ( |2.7p for F 2 ns- The left part is analogous to Fig. 10, in the right part the resulting 
a 5 s corrections to the result in Fig. 4 are compared to the a\ (N 3 LO) contribution. The 
scales of the graphs are the same as in Fig. 9. 
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represent good approximations at both orders. The Pade approximations are somewhat 
smaller, however, this offset seems to decrease with the order in a s . In the left part of 
Fig. 11 we present the corresponding N 4 LO predictions. The inner three curves have been 
derived from the central N 3 LO results of section 5. The PMS and ECH are again very 
similar, they are not shown separately. The impact of the present uncertainty of the N 3 LO 
kernels, dominated by the estimate ( ^3|) of the four-loop splitting functions, is included 
in the two dotted curves which represent our final estimate for the N 4 LO term and its 
uncertainty. In the right part of Fig. 10 the N 4 LO corrections to the results of Fig. 4 
are compared to the N 3 LO contribution. Within the large uncertainties of the latter, 
these results are consistent with the NNLL soft-gluon prediction shown in Fig. 9. The 
consequences of including the N 4 LO estimates have been presented in Fig. 7 and Eq. ( |5.5|) , 
respectively, for the renormalization-scale stability and the determination of a s . 

Finally the «' >5 infinite-order predictions of the soft-gluon resummation (|2.17|) (using 
the minimal prescription contour |27|) and the Pade approximations ( |6.4[ ) are compared 
in Fig. 12. For the present uncertainties (|3.14| ) — the curves in the figure refer to £ = 8 



(upper), D 2 = in Eq. (|2.18|) (middle) and £ = 12 (lower) — it is not possible to draw any 



conclusions from the soft-gluon result. The Pade summation, on the other hand, provides 
rather definite predictions: The terms beyond can be expected to have a very small 
impact at x < 0.75. For our standard reference value a s = 0.2 their effect reaches about 
the size of the N 4 LO and N 3 LO contributions at x ~ 0.9 and x ~ 0.95, respectively. 
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Figure 12: Predictions of the soft-gluon resummation (|2.17|) and the Pade approximations 



for the contributions beyond N LO to the logarithmic scale derivative of F 2 = F<< 
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7 Summary 



We have investigated the predictions of massless perturbative QCD for the scaling viola- 
tions of the most important non-singlet structure functions in unpolarized DIS, extending 
our previous NNLO results |L1] to N 3 LO and N 4 LO for the region x > 10 -2 . The main 



objective of this extension is to reduce the theoretical uncertainty of determinations of 
a s from inclusive DIS to about 1%, an accuracy which is sufficient to make full use of 
present and future structure function measurements. Our results also facilitate improved 
determinations of power-suppressed contributions to the structure-function evolution by 
fits to data, especially at large x where the uncertainties are still sizeable at NNLO. 

The new ingredients entering the N 3 LO physical evolution kernels are the four-loop 
splitting functions and the three-loop coefficient functions. The impact of the former 
quantities is expected to be very small in the MS scheme at x> 10~ 2 ; it has been estimated 
by a Pade approximation assigned a 100% uncertainty. For the latter quantities we 
have derived approximate expressions based on the available integer moment |7], |S| and 



soft gluon [15|, [Tc]] results. The effect of these functions is very well under control at 
x > 10 -2 , and almost perfectly at x > 0.1. In fact, the uncertainty of the splitting functions 
dominates the small residual uncertainties of the evolution kernels. Hence the accuracy 
of our present N 3 LO results will be superseded only by a future four-loop calculation. 



We have also studied the predictions of the NNLL soft-gluon resummation |T7[ and 
of the Pade, PMS and ECH approximations [|TH], |T!], [2"D]| . Presently the predictions of the 
resummation for the physical evolution kernels beyond N 3 LO (in any case applicable only 
at x > 0.8) suffer from the incomplete determination of the soft-gluon parameters B 2 and 
D2, a problem which will be removed by forthcoming exact calculation of the three- loop 



coefficient functions [Q. The Pade, PMS and ECH approximations are found to agree 
rather well with the NNLO and N 3 LO results for the evolution kernels; these approaches 
seem to provide reliable predictions of the effects at N 4 LO and beyond. 

For a s < 0.2 the N 3 LO and N 4 LO corrections at /i r ~ Q are very small at x<0.6 and 
x < 0.8, respectively, especially for the most accurately measured structure function F 2 . 
Consequently the central values of a s determined from the non-singlet scaling violations 
hardly change any more once the larger NNLO terms have been included, see also ref . @ . 
The scale uncertainty of the resulting o; s (M|) is reduced to the unproblematic level of less 
than 1% at N 3 LO and 0.5% at N 4 LO. In order to ensure an overall theoretical accuracy 
of about 1% also the heavy quark (especially charm) mass effects need to be controlled 
with this precision. We will address this point in a forthcoming publication. 

Fortran subroutines of our approximations of the three-loop coefficient functions in 
section 3 and of the parametrizations of the convolutions entering the evolution kernels 
in section 4 can be found at [http : //www . lorentz . leidenuniv . nl/^avogt . 
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Appendix: Convolutions of +- distribution 

The convolutions T>i ® T>j of the -(—distributions ( |3.2|) for i + j < 4 are given by 
V ®V = 2V 1 -( 2 d(l-x) 
V X ®V = -V 2 -( 2 V + ( 3 d(l-x) 

V X ®V X = V 3 -2( 2 V 1 + 2( 3 V -^( 2 2 6(l-x) 

V 2 ®V = -V z -2$ 2 V 1 + 2&Vv-\(*8(l-x) 
6 5 

V 2 ®V X = ^V 4 -3( 2 V 2 + Q( 3 V 1 -( 2 2 V + (4( 5 -2( 2 ( 3 )5(l-x) 
V 2 ®V 2 = ^P 5 -4C 2 P 3 + 12C3^2-4C 2 2 ^i + (l6C5-8C 2 C3)^o (A.l) 
+ (4C 3 2 - ^ Cl) 5(1 -x) 

V 3 ®V = ^-V 4 -3( 2 V 2 + 6( 3 V 1 -^( 2 2 V + 6( 5 5(l-x) 
4 5 

o 97 
V 3 ®V 1 = -V 5 -4( 2 V 3 + 12( 3 V 2 - _C 2 2 P 1 + (l8C5-6C 2 C 3 )^o 

+ (3C|-|C 2 3 ) 5(1 -x) 

6 48 192 

V 4 ®V = -V 5 -4( 2 V 3 + 12( 3 V 2 - —(IV 1 + 24( 5 V -—(l5(l-x) 

up to integrable contributions dealt with numerically in section 4. Here Q stand for the 
Riemann ^-function, and (4 and (q have been expressed in terms of ( 2 and £3, respectively. 
A convenient method to derive (and to extend, if required) eqs. (|A.1| ) is by using the 
relation between T>i and the harmonic sums discussed, for example, in ref. |34] . 
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